Scaling and correlations in the dynamics of forest-fire occurrence 
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Forest-fire waiting times, defined as the time between successive events above a certain size in a given region, 
are calculated for Italy. The probability densities of the waiting times are found to verify a scaling law, despite 
that fact that the distribution of fire sizes is not a power law. The meaning of such behavior in terms of the 
possible self-similarity of the process in a nonstationary system is discussed. We find that the scaling law arises 
as a consequence of the stationarity of fire sizes and the existence of a non-trivial "instantaneous" scaling law, 
sustained by the correlations of the process. 



In the last years, many natural hazards, like earthquakes, 
volcanic eruptions, landslides, rainfall, solar flares, etc., 
and other similar-in-spirit phenomena in condensed-matter 
physics have been shown to be characterized by a power-law 
distribution of event sizes, over many orders of magnitude in 
some cases [????]. This kind of distribution has profound 
implications for the nature of these phenomena, as it indicates 
that extreme events do not constitute a case separated from 
the smaller, ordinary ones; rather, the events are generated by 
a mechanism that operates in the same way for all the dif- 
ferent scales involved, and a characteristic size of the events 
cannot be defined. In this way, a reasonable question such as 
"which is the typical size of the earthquakes in this region?" 
is impossible to answer. 

Comparison with simple self-organized-critical (SOC) 
cellular-automaton models suggests that the events that define 
the dynamics in these phenomena consist of a small instabil- 
ity or excitation that propagates as a very rapid chain reaction 
or avalanche through a medium that is in a very particular 
state, similar to the critical points found at continuous phase 
transitions in condensed-matter physics [????]. The dis- 
sipation produced by each avalanche would act as a feedback 
mechanism that balances a slow energy input and maintains 
the system close to the critical state. 

Of special interest is the case of forest fires, for which 
cellular-automaton models yielded a power-law behavior for 
the distributions of burned areas (which are a measure of the 
size of the events), and showed the previous mechanisms at 
work [? ? ]; curiously, it was not until much later that Mala- 
mud et al. observed power law distributions for real forest 
fires, with exponents around 1.4 for the probability density 
(i.e., non-cumulative distribution) of fire sizes [???]. 

Nevertheless, this issue is still open, as other studies with 
different data do not agree with a simple power-law behav- 
ior: Ricotta et al. [? ] postulated that fires of large sizes, 
due to negative economic and social effects, are reduced by 
the massive human intervention; therefore, less than expected 
large fires occur, leading to an increase (in absolute value) of 
the power-law exponent in that regime. Reed and McKelvey 
[? ], using the concept of extinguishments growth rate, pre- 
sented a four-parameter "competing hazards" model provid- 
ing the overall best fit. In a subsequent paper, Ricotta et al. [? 
] have observed that a multiple power-law behavior, denoted 



by the presence of different power-law ranges delimited by 
cut-offs, is due to dynamical changes, linked to "more or less 
abrupt changes in the landscape-specific process-pattern in- 
teractions that control wildfire propagation, rather than statis- 
tical inaccuracies". Therefore, the appearance of different size 
ranges with different power-law exponents can be accounted 
for different dynamics, involving topographic, climatic, vege- 
tational, and human factors [? ]. 

In addition to the size of the events, the dynamics of event 
occurrence is of fundamental interest. Notably, the tempo- 
ral properties of some popular SOC cellular-automaton mod- 
els were shown to be described by a trivial Poisson process, 
which prevented progress in this aspect until very recently, 
when it has been concluded that this picture is not appropri- 
ate [? ]; in parallel, it has been found that real systems show 
a very rich behavior in time, with power-law distributions for 
the time between events and/or scaling laws for these distribu- 
tions [???]. The existence of such scaling laws has implica- 
tions no less deep than the fact of having a power-law distribu- 
tion of event sizes, although they have been much less studied: 
(i) the scaling law reflects the fact that the occurrence of large 
events mimics the process of occurrence of smaller ones (and 
this behavior is not implicit in the distribution of event sizes), 
thus allowing to model the scarce big events on the basis of the 
abundant small ones; (ii) the scaling law is the signature of the 
invariance of the process under a renormalization-group trans- 
formation, which strengths the links between natural hazards 
and critical phenomena [? ]. 

We study in this paper the relation between the temporal 
properties of forest-fire occurrence and the size of the fires, 
using the AIB (Archivio Incendi Boschivi) fire catalog com- 
piled by the Italian CFS (Corpo Forestale dello Stato) for all 
Italy [? ], covering the years 1998—2002 (included) and con- 
taining 36821 fires. In order to characterize the overall behav- 
ior, we measure for the whole catalog the probability density 
of the burned areas s, defined as 



Prob [s < area < s + ds] 

D(s) = 

as 



(1) 



where ds is the bin size (small enough to sample almost con- 
tinuously D(s) but large enough to guarantee statistical signif- 
icance); the resulting shape for D(s) is shown in Fig. [T] Al- 
though a power law could be fit to the data, it is clearly seen 
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that the curve is continuously bending downwards, which is 
the characteristic of a lognormal distribution, 

(lns-/i) 2 \ fe»\ l+ ^- 

(2) 

with jj. and a the mean and standard deviation of In s, and C 
a correction to normalization due to the fact that the fit is not 
valid for all s. In this way, for each 2a 2 that In s is away from 
ji the exponent of the previous pseudo-power law increases in 
one unit (in other words, each decade s is above e' 1 increases 
the exponent in In 10/ (2a 2 )). When s is measured in hectares 
(ha), the results of the best fit yield ji = —0.35, a = 9.5, and 
C — 6.7; this fit holds not only for the full data but it can be 
verified that also describes smaller parts of the country and 
shorter periods of time. In any case, we have no means to 
conclude if the deviation from a power-law behavior is due to 
human extinction efforts or to the territorial characteristics of 
a high-populated country. 

From the distribution of sizes, knowing the total number 
of events, it is possible to calculate the mean waiting time 
(or recurrence time) for events above a certain size s c [? ]; 
however, looking at the individual values of the waiting times 
one sees that they are broadly distributed and therefore the 
mean values are not very informative about the dynamics of 
the process; so, in order to investigate the temporal properties 
of fire occurrence it is necessary to look at the whole waiting- 
time distribution. To be precise, the procedure is as follows: 
once a spatial area, a time period, and a minimum event size, 
s c , are selected, the fire history is described as a simple point 
process, {?o,fi,f2 ■ ■ where t\ denotes the time of occurrence 
of fire i. For this process, the set of waiting times, defined 
as the time intervals between consecutive events, is obtained 
straightforwardly as t, = f, — i . Important insight into the 
nature of the process may be obtained by considering s c not 
as a constant but as a variable parameter [? ? ], and then, 
the waiting-time probability density for the selected window, 
defined in the same way as in Eq. (fTJ, will be also considered 
as a function D(t;s c ) of the minimum size s c . 

For the whole country and the total temporal extension of 
the catalog we obtain the different set of curves displayed 
in Fig. 0a). We might fit a (decreasing) power law for 
each distribution, but the exponent would decrease with the 
increase of the minimum size s c . Instead, it is more con- 
venient to rescale the distributions in order that all of them 
have the same mean and can be properly compared; this is 
accomplished by the scale transformation T — > R(s c )t and 
D(x;s c ) — > D(t,s c )/R(s c ), where R(s c ) is the rate of fire oc- 
currence, defined as the mean number of fires per unit time 
with s > s c (that is, the inverse of the mean of each distribu- 
tion). The results of the rescaling, as shown in Fig. |2jb), lead 
to a collapse of the rescaled distributions into a single function 
F, signaling the fulfillment of a scaling law, 

D(t;sc)=R(s c )F(R(s c )t), (3) 

in the same way as for several natural hazards [? ? ? ] and 



other avalanche-like processes [? ? ]. 

The rescaled plot unveils more clearly the behavior of the 
distributions: instead of different power laws, what we have 
is a unique shape, but at different scales. Again, the apparent 
continuous decrease of the exponent with the rescaled time, 
9 = R(s c )x, suggest a lognormal shape for F(9) as that of Eq. 
(fJJ, where now we will use tildes to denote the parameters. 
The best fit yields ft, = — 2.0 and a = 2.0, fixing C = 1 . Notice 
that now we have the constraint that the mean of the rescaled 
distribution, 0, has to be one; as 6 = e^+ a I 2 , this leads to 
PL = -a 2 j2. 

It is remarkable that, unlike earthquakes, solar flares, or 
fractures [? ? ? ? ], forest fires fulfill a scaling law 
for the waiting time distributions without displaying power- 
law distribution of event sizes. We could conclude that we 
have self-similarity in size-time without having scale invari- 
ance in size alone. This self-similarity means that for the 
linear scale transformation T — > ax and s c — > bs c , the value 
of b which guarantees scale invariance is given by R(bs c )a = 
R(s c ), which means that b does not only depend on a, as in the 
case of a power-law distribution of sizes, but it also depends 
on s c . This would be equivalent to define an artificial new size 
variable enforcing that it be power law distributed. However, 
although this picture describes a kind of self-similarity, it is 
not a sufficient condition. Indeed, the seasonality of fire oc- 
currence prevents self-similarity in size-time: five years of fire 
occurrence cannot be equivalent to one year of smaller fires, 
as there is a clear annual modulation in fire occurrence; nev- 
ertheless, for a fixed time window still the small events are a 
model for the occurrence of the big ones. 

Which is then the origin of the scaling law (0? It is not 
difficult to relate it with the stationarity of fire sizes and with 
the existence of a scaling law for the "instantaneous" waiting- 
time distributions. Indeed, D(t;s c ) is a statistical mixture 
of those instantaneous waiting-time distributions D t (x;s c ), 
which, when the scale of variations of the rate is much larger 
than the corresponding mean waiting time, take into account 
that fire occurrence is not stationary but change with time 
f; if it is only the instantaneous rate r(t;s c ) (defined as the 
number of fires per unit time in a small time interval around 
t) what determines fire occurrence, we can write D t (x;s c ) = 
D(v,s c \r(t;s c )) and then, 

I [-''max 

D(x;s c ) = — — - / rD(T;s c \r)p(r;s c )dr, 

R ( s c) Jr mm 

where p(r;s c ) is the density of rates, i.e., the fraction of the 
time the rate is in a particular small range of values, divided 
by that range [? ]. Assuming the stationary nature of fire 
sizes (notice that this is not incompatible with the nonstation- 
arity of time occurrence), this means that r(t;s c ) = pr(t;so), 
where the fraction p is the probability of having a fire larger 
than s c knowing that it has been larger than so, p = Prob[s > 
s c ]/Prob[i > So]; this implies that the density of rates fulfills a 
scaling law, p(r;s c ) — p~ 1 p(p~ l r;so) = p~ l g(p~ l r). Finally, 
with the hypothesis that D(x;s c \r) verifies as well a (instanta- 
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neous) scaling law, D(x;s c \r) — rf(rx), we get 
1 fpb 

D(r,s c ) = — rf(rT)p- l g(p- l r)dr, 

pl<0 J pa 

with r mjn = pa, r max = pb, and R(s c ) = pR(s ) = pRo- A sim- 
ple change of variables reveals that D(x;s c ) is a function of 
the form pF(px) = p x 2 f(pxx)g(x)dx, which is equivalent 
to the scaling law (0. In other words, if fire occurrence under 
hypothetical stationary conditions verifies a scaling law for the 
waiting times (which in this case would be a reflection of the 
self-similarity of the stationary process, as explained above), 
non-stationary conditions keep that scaling valid (with a dif- 
ferent scaling function) as long as fire size remains stationary 
and the rate does not become too small for this description to 
be invalid. [On the other hand, for rates so small that the mean 
waiting time is much larger than the larger scale of variation 
of the rate itself (whose existence is not known), the structure 
of r(t) would become irrelevant and the waiting-time distri- 
bution would tend to the exponential form characteristic of 
Poisson processes.] 

In order to support our argument for the existence of the 
scaling law (01 we show in Fig. [3]the stationarity of fire sizes, 
by means of the evolution of r{t\s c ) for different s c , and how 
the different curves collapse when they are rescaled by their 
mean, R(s c ); it is also easy to check that the distribution of 
rates verifies a scaling law. The last hypothesis, the scaling 
of D(x;s c \r) is more difficult to demonstrate due to the daily 
oscillations of r(t;s c ), which makes that the rate can be con- 
sidered approximately constant only for a few hours, corre- 
sponding to those of the daily maximum and minimum hazard 
(between 1 p.m. and 4 p.m. and between 1 a.m. and 9 a.m., 
respectively). This short range of variation leads to very low 
statistics; nevertheless, for the periods of the year of maxi- 
mum fire occurrence (for about one month in the summer) the 
maximum and minimum daily rates are fairly constant for dif- 
ferent days, which allows to improve the statistics. The results 
obtained in this way are shown in Fig. |U although they are not 
conclusive. Essentially, they are compatible with an instanta- 
neous scaling law, with perhaps an exponential instantaneous 
distribution, D(x;s c \r) ~ re~' T , but the statistical errors are 
large; in any case, the hypothesis of the instantaneous scaling 
law cannot be rejected. 

If we find an exponential form for the instantaneous distri- 
butions, does this mean that the dynamics can be described by 
a nonstationary Poisson process? This is the simplest model 
for nonstationary behavior, for which the events take place at a 
rate that does not depend on the occurrence of the other events, 
as in the simple (stationary) Poisson process, but with the dif- 
ference that the rate changes with time (independently on the 
process, we can imagine the rate is related to the meteorolog- 
ical conditions, not affected by the presence of fire or not). 
This leads indeed to exponential instantaneous distributions 
(provided the rate is not too small), although the reciprocal is 
not true, in general. If, in addition, the size of the events con- 
stitutes an independent random process, this ensures the exis- 
tence of a scaling law for the instantaneous distributions (as 



Poisson processes are invariant under random thinning plus 
rescaling, see [? ]). The nonstationary Poisson process has 
been recently used for earthquake occurrence, see Ref. [? ]. 

A test to verify if a process is of the nonstationary Poisson 
type was introduced by Bi et al. [? ]. One only needs to 
compute for each i the statistics hj = 2T m , M ;/(2T m ,„/ + T n eigi), 
where T m / ni is the minimum of T, and T,- + i, and T ne i g i is the 
length of the interval neighbor of the minimum one opposite to 
the one used in the comparison, i.e., T,-_i or T, + 2 respectively. 
Under the hypothesis we want to test, both X ne i g and 2T m ;„ are 
independent and exponentially distributed with approximately 
the same rate, r(f), and therefore it can be shown that h is 
uniformly distributed between and 1 . 

The application of the test to the fire data yields catastrophic 
results, see Fig. [5] The obtained probability density for h is far 
from uniform, with very large peaks for precise /z-values. This 
is due to the discretization of fire occurrences in the catalog, 
which are determined verbally and therefore rounded mainly 
in units of 10 or 15 min; this favors particular values of x and 
there fore of h (2/3, 4/5, 1/2, etc.). We can correct this effect 
by the addition of a uniform random value between -5 min and 
5 min to each occurrence time U, then the peaks in the distribu- 
tion of h disappear and its shape gets closer to a uniform one; 
however, the difference is significant. We have verified that 
the difference is not due to the random addition we have per- 
formed: simulation of a non-stationary Poisson process where 
the occurrences are rounded in intervals of 10 min yields a per- 
fect uniform distribution when this discretization is corrected 
by the uniform random addition just explained (Fig. [5]). In 
consequence, this model does not seem suitable for fire occur- 
rence, and although the instantaneous distributions are close to 
exponential (Fig. HI, this is not a sufficient condition to have a 
nonstationary Poisson process, as the absence of correlations 
is equally important for it. 

If we reject the nonstationary Poisson process with inde- 
pendent sizes as a model of fire occurrence, the only way to 
get a scaling law for the instantaneous process is by means of 
orchestrated correlations between sizes and occurrence times 
[? ]. In order to establish the existence of such correlations 
we proceed to study conditional size distributions, defined as 
in Eq. (Q~|l but with an additional condition for the compu- 
tation of the probability. We consider D(s\s„ re > s[.), which 
accounts for the size of the events for which the size of the im- 
mediate previous-in-time event, s pre , is above a given thresh- 
old s' c . The results in Fig. |6]show that an increase of s' c triggers 
a greater proportion of large fires, i.e., large fires are followed 
by large fires. The dependence of a fire size on the previous 
size is small but significant, unlike to what happens for earth- 
quakes, where correlation between their magnitudes has not 
been detected [? ? ] (nevertheless, for an alternative view see 
Ref. [? ]). On the other hand, the dependence of waiting times 
on the size of the event defining the starting of the waiting pe- 
riod can be measured by D(x;s c | s pre > s' c ) showing how large 
fires cause a decrease in the number of long recurrence times, 
i.e., those fires tend to be closer in time to the next fires. The 
effect is again small, but clearly detectable, and in this case 
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has a counterpart for earthquakes, where the Omori law for 
aftershocks implies the same behavior. 

But the correlations between fires are not only with the pre- 
vious event; its range can be quantified by means of the fol- 
lowing auto-correlation function, 

c(j;s c ) = ((logs; - £)(logs i+j ~ I)) al 2 , 

where I is the arithmetic mean of the logarithm of the size 
(i.e., the logarithm of the geometric mean of the size), and 
Ot> is the standard deviation of the logarithm; both i and Oi 
depend on s c . Notice that although the process is not station- 
ary, the stationarity of the size gives sense to the autocorrela- 
tion function defined in this way. The results for this function 
are shown in Fig. UJ and compared with the same correlation 
function calculated for a reshuffled version of the catalog, for 
which the size of the events are randomly permuted, breaking 
the correlations between them (which should yield an autocor- 
relation function fluctuation around zero). The conclusion is 
that positive correlations extend significantly beyond several 
hundreds of events (for events of size larger than 1 ha). 



More clear is the behavior of the autocorrelation as a func- 
tion of time; as the process is not stationary both functions are 
not equivalent. We define 

c(A;s c ) = ((logs(t)-i)(log S (t+A)-I))oi 2 , 

where s(t ) denotes the size of the fire that happens at time t 
(we slightly change notation, for convenience). The average 
is taken over all times t and t + A for which there are fires, 
this yields the results of Fig. UJ The correlation is again posi- 
tive, but larger in this case, suggesting that real time is a better 
variable to describe the evolution of correlations, which ex- 
tend for about 10000 min, i.e., roughly 1 week. It is likely 
that these correlations are mediated through the meteorologi- 
cal conditions. 

In summary, the dynamics of forest-fire occurrence shows 
a complex scale-invariant structure at any time, modulated by 
seasonal and daily variations and orchestrated by means of 
broad-range correlations. 
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FIG. 1: Probability density of fire sizes in Italy, from 1998 to 2002 (included). The error bars are calculated for one standard deviation in the 
number of counts. The fit is the lognormal distribution whose parameters are given in the text. 




r (min) R{s c )t (dimensionless) 



FIG. 2: (a) Probability densities of fire waiting times for Italy, from 1998 to 2002 (included) for different minimum burned areas (s c = 1 ha 
to s c = 300 ha); waiting times smaller than 5 min are not plotted. The slope of the tail seems to decrease with increasing s c . (b) The previous 
densities after rescaling by the mean fire rate, adding error bars (corresponding to one standard deviation). The data collapse indicates the 
existence of a scaling law (see text) and allows a unified description of the shape of the density, in terms of a lognormal distribution, rather 
than as a power law. 
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FIG. 3: Monthly rates of forest-fire occurrence in Italy as a function of time, for events larger or equal than s c , with s c ranging from 1 ha to 30 
ha (top curves). In addition, the rescaling of the rate by their mean R(s c ) is shown (bottom collapse of curves), indicating the stationarity of 
fire sizes as well as the nonstationary occurrence in time. 




FIG. 4: Rescaled distributions of waiting times during annual maximum rate periods, corresponding to (in years) 1998.50-1998.65, 1999.60- 
1999.66, 2000.56-2000.67, and 2001.54-2001.64 (same symbols) separated for daily minimum rate, 1 a.m. - 9 a.m. and maximum, 1 p.m. - 
4 p.m. Different minimum sizes are used. The exponential scaling function is shown for comparison. 
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FIG. 5: Probability densities of the statistics h, for the original catalog, for the smoothed catalog (corrected), and for a simulation of a 
nonstationary Poisson process with the same rate as the real process to which the same discretizaion and smoothing procedure has been 
applied. 
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FIG. 6: Conditional probability densities of fire sizes in Italy when the previous fire is larger than s' c . Only fires greater than 1 ha have 
been taken into account. The case s pre > 1 ha corresponds to the unconditioned distribution shown in Fig. 1 . The change of the conditional 
distributions with respect the unconditioned case is an indication of positive correlations between fire sizes. 




FIG. 7: Autocorrelation functions for the Italian catalog and for a version with reshuffled sizes. In the latter case it is the modulus of the 
autocorrelation what is shown, as the function fluctuates around zero and it is equally likely that it is positive or negative. 



